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SUMMARY 

A  numerical  method  is  developed  for  the  solution  of  Stokes’s  stream 
function.  The  method  is  applied  to  the  existing  Low  Speed  Wind  Tunnel 
(LSWT)  contraction  shape.  The  calculated  wall  pressure  coefficient 
distribution  is  compared  to  experimental  measurements  to  test  the 
approximation  of  axial  symmetry  made  for  the  octagonal  sections  and  for  the 
effects  of  neglecting  viscous  forces.  Reasonable  agreement  is  obtained.  The 
existing  contraction  is  known  to  suffer  from  intermittent  boundary  layer 
separation  near  the  entry  caused  by  adverse  pressure  gradients.  A  new 
contraction  shape  is  suggested  that  uses  the  contraction  ratio  of  the  existing 
LSWT  so  that  the  new  design  can  be  installed  within  the  existing  LSWT  in  the 
form  of  an  inner  skin.  The  new  design,  which  is  about  10%  longer  than  the 
existing  contraction,  has  an  adverse  pressure  gradient  near  the  inlet  which  is 
about  1/5  of  that  of  the  existing  contraction  while  the  adverse  pressure 
gradient  near  the  exit  is  around  the  same  value.  It  is  expected  that  the 
reduced  adverse  pressure  gradient  near  the  entry  will  be  small  enough  to  allow 
the  inlet  boundary  layer  to  remain  attached  thus  leading  to  a  substantial 
improvement  to  the  working  section  flow  quality.  The  new  design  could  be 
installed  upstream  of  the  existing  contraction  to  provide  increased  working 
section  length. 
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1.  INTRODUCTION 


The  acceleration  'of  the  flow  in  a  wind  tunnel  contraction  prior  to  the 
working  section  ensures  that  the  rest  of  the  tunnel  circuit  is  subject  to  reduced 
losses  in  the  tunnel  wall  boundary  layers,  turning  vanes,  screens  and  honeycomb. 
Mean-flow  non-uniformities  are  reduced  so  that  a  more  uniform  velocity  profile  is 
obtained  in  the  working  section.  The  streamwise  turbulence  intensity  in  the 
working  section  is  also  reduced  by  stretching  of  the  wakes  from  the  turning  vanes, 
honeycomb  and  screens. 

Since  the  flow  is  accelerated  by  the  contraction  the  pressure  gradients  on 
the  walls  are  generally  favourable.  However  it  can  be  shown  that  for  a  finite 
length  contraction  two  regions  of  adverse  pressure  gradient  exist  on  the  walls  in 
the  vicinity  of  the  entrance  and  the  exit.  The  presence  of  the  contraction  causes 
curvature  of  the  streamlines  in  the  upstream  section  of  constant  cross-sectional 
area  causing  the  velocity  to  fall  near  the  wall.  The  wall  pressure  continues  to 
rise  and  reaches  a  local  maximum  near  the  inlet  (but  within  the  contraction) 
before  beginning  to  fall.  Similarly  a  pressure  minimum  (undershoot)  occurs  on  the 
wall  near  the  exit  of  the  contraction.  As  a  result  the  velocity  profile  across  the 
exit  plane  of  a  contraction  is  non-uniform  and  in  a  short  contraction  the 
non-uniformity  may  be  as  high  as  10%.  The  wall  velocities  approach  the  mean 
value  asymptotically  so  that  a  model  under  test  in  the  working  section  should  be 
placed  some  distance  downstream  from  the  contraction  exit,  i.e.  the  useful  length 
of  the  working  section  is  reduced.  The  two  regions  of  adverse  pressure  gradient 
can  have  a  further  detrimental  effect  by  causing  large-scale  intermittent 
boundary  layer  separation  resulting  in  increased  unsteadiness  of  the  flow  in  the 
working  section. 

The  ARL  low  speed  wind  tunnel  (LSWT)  has  been  in  almost  continuous 
operation  since  it  was  commissioned  in  1941.  Figure  1  is  a  sketch  of  the  tunnel 
showing  the  important  dimensions.  The  contraction  is  characterised  by  an 
initially  rapid  rate  of  contraction  followed  by  a  much  more  gradual  variation  near 
the  exit.  It  appears  that  the  designers  were  more  concerned  about  boundary  layer 
separation  caused  by  adverse  pressure  gradients  near  the  contraction  exit  rather 
than  near  the  entry.  Batchelor  and  Shaw  U944>  found  that  the  total  head  in  the 
working  section  fluctuates  by  around  1%  throughout  the  test  section.  Larger 
fluctuations  of  up  to  10%  were  experienced  near  a  fillet.  Turbulence  intensity 
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measurements  by  D.C.Collis  and  M.J.  Williams  in  1946  indicate  that  the 
streamwise  turbulence  intensity  varies  from  0.2  to  0.7%  depending  on  the  location 
in  the  working  section.  The  numerical  computations  of  Batchelor  and  Shaw 
indicate  that  a  major  source  of  the  unsteadiness  may  be  due  to  intermittent 
boundary  layer  separation  near  the  entry  of  the  contraction  caused  by  adverse 
pressure  gradients.  It  is  now  well  known  (see  Bradshaw  1973)  that  boundary 
layers  are  more  unstable  over  concave  surfaces  (e.g.  near  a  contraction  inlet), 
than  those  over  convex  surfaces  (e.g.  near  a  contraction  exit).  Recent  wool  tuft 
flow  visualisation  studies  have  confirmed  that  the  boundary  layers  at  the  entry  to 
the  contraction  suffer  from  intermittent  separation.  The  present  configuration 
of  the  LSWT  also  has  an  unusually  short  working  section  which  has  required  testing 
to  be  performed  in  regions  that  may  be  too  close  to  the  contraction  exit  for 
adequate  flow  uniformity. 

Until  recently  theoretical  methods  of  contraction  design  were  considered 
unsatisfactory  and  most  designers  relied  heavily  on  empirical  methods.  The 
analytical  methods  produce  series  solutions  which  are  of  little  practical  use. 
Some  of  the  series  solutions  are  inappropriate  since  they  are  periodic  along  the 
contraction  axis.  Others  require  impractically  long  axial  lengths  to  obtain 
satisfactory  entry  and  exit  conditions.  More  recently  a  number  of  numerical 
solutions  have  become  available  for  guiding  the  designer.  The  advantage  of  using 
numerical  methods  is  that  solutions  may  be  obtained  for  an  arbitrary  contraction 
shape.  Recent  studies  have  been  made  for  families  of  contraction  shapes  that  are 
specified  by  a  simply  polynomial  relating  the  radius  to  axial  distance  (Morel 
1975).  However  there  has  been  little  published  experimental  evidence  to  confirm 
the  adequacy  of  these  designs.  Willis  and  Hunt  (1982)  suggested  a  new  design  for 
the  LSWT  where  the  contraction  ratio  was  increased  from  4:1  to  7:1. 
Implementation  of  their  design  would  require  major  modifications  to  the  steel 
shell  of  the  I-SWT  which  would  be  an  expensive  and  time  consuming  task.  In  this 
report  a  numerical  method  is  developed  for  solving  Stokes’s  stream  function.  On 
the  basis  of  results  from  the  calculation  method  a  new  LSWT  contraction  design  is 
suggested.  The  contraction  ratio  is  4:1  so  that  the  new  design  could  be  installed 
within  the  existing  LSWT  steel  shell. 
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2.  GOVERNING  DIFFERENTIAL.  EQUATION 

The  continuity  equation  in  an  incompressible  axisym  metric  flow  is  given  by, 


i  3(rV)  3U  _  0 
r  ar  +  ax 


where  x  and  r  are  coordinates  in  the  axial  and  radial  directions  and  U  and  V  are 
the  respective  fluid  velocities. 


This  equation  is  satisfied  by  defining  a  stream  function  0  such  that 


U  =  j;  and 
r  3r 


_L  Ml 
r  ax 


For  irrotational  flow,  —  -  —  =  0 

so  that,  £*  +  0  (1) 

ax2  r  ar  ar2 

The  stream  function  ^satisfying  this  equation  is  known  as  Stokes’s  stream 
function.  Suitable  boundary  conditions  for  contraction  are  : 

(1)  0=0  along  the  axis. 

(2)  0  =  =  constant,  along  the  boundary. 

2 

(3)  0 =  (r/r  )  0  across  the  entry  to  the  region  of  constant  cross-sectional 

area  upstream  of  the  contraction  (rw  =  radius  of  tunnel 

2 

boundary).  This  gives  U  =  20^/r^  =  constant. 

Provided  that  the  length  of  the  parallel-walled  entry 
region  is  large  enough  then  the  radial  velocity  component 
will  also  be  zero  i.e.  the  entry  flow  will  be  uniform  and 
parallel  to  the  axis. 
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(4 )  3  ip/3x  =  0  across  the  exit  of  the  region  of  constant  cross-sectional  area 

downstream  of  the  contraction. 

3.  THEORETICAL  METHODS  OF  CONTRACTION  DESIGN 

Most  analytical  methods  for  contraction  design  are  considered  inappropriate 
for  practical  use.  By  far  the  most  common  analytical  technique  has  been  the 
method  of  separation  of  variables  of  equation  1.  Depending  on  the  choice  of 
separation  constant  three  possible  series  solutions  may  be  constructed.  Each  of 
these  solutions  suffers  from  serious  shortcomings.  The  first  type  of  solution  is 
restricted  to  families  of  axial  velocity  distributions  which  permit  analytical 
expression  of  their  derivatives  e.g.  Tsien  (1943),  Szceniowski  (1943)  and  more 
recently  Cohen  and  Ritchie  (1962).  The  major  difficulty  in  applying  this  type  of 
solution  arises  from  the  stringent  requirements  on  the  form  of  the  axial  velocity 
distribution.  Barger  and  Bowen  (1972)  generalised  the  method  to  cater  for  a 
wider  range  of  velocity  distributions  so  that  the  method  could  be  applied  to  a 
greater  variety  of  duct  shapes.  Nevertheless  their  approach  is  restrictive  in  the 
sense  that  only  a  limited  number  of  contraction  shapes  can  be  obtained  from  the 
solutions.  The  velocities  at  the  entrance  and  the  exit  are  neither  uniform  nor 
parallel  and  a  rather  long  contraction  is  required  to  avoid  a  strong  adverse 
pressure  gradient  at  the  entrance.  The  second  type  of  solution  is  periodic  in  the 
streamwise  direction.  Thwaites  (1946)  obtained  series  solutions  in  terms  of 
Bessel  functions.  Since  eqation  (1)  is  elliptic  the  effect  of  periodic  repetition  of 
the  contraction  shape  will  produce  a  flow  different  to  that  obtained  with  upstream 
and  downstream  regions  of  constant  area.  It  is  possible  to  specify  parallel 
entrance  and  exit  flow  but  the  velocity  distributions  will  not  be  uniform. 
Alternatively  a  uniform  velocity  distribution  may  be  specified  but  at  the  expense 
of  parallel  entry  and  exit  flow.  However  it  can  be  shown  in  general  that  the 
entrance  and  exit  velocity  distributions  cannot  be  both  uniform  and  parallel  if  the 
entrance  and  exit  are  a  finite  distance  apart.  The  third  type  of  solution  is  in  the 
form  of  an  a  periodic  (exponential)  series  which  produces  solutions  for  the 
contraction  shape  which  are  infinitely  long.  A  practical  design  based  on  this  type 
of  solution  would  need  to  be  truncated  which  may  produce  unknown  and 
undesirable  affects.  Other  techniques  that  have  been  tried  include  Whitehead  and 
Wu  (1951)  who  used  an  approximate  method  where  the  velocity  potential  and 
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stream  function  of  a  two  dimensional  contraction  are  used  as  independent 
variables  and  the  equation  for  the  Stokes’s  stream  function  is  solved 
numerically.  Smith  and  Wang  (1944)  used  a  distribution  of  singularities  such  as 
ring  vortices,  sources  and  sinks  to  design  a  family  of  contra ctions. 

It  is  unfortunate  that  there  has  been  little  experimental  confirmation  of  the 
adequacy  of  these  analytical  solutions.  One  exception  is  the  work  by  Bossell 
(1969)  who  built  a  16:1  contraction  based  on  a  periodic  solution  of  Stoke’s  stream 
function.  Unfortunately  his  apparatus  was  rather  small  (exit  diameter  of  150mm) 
so  that  the  experimental  velocity  profiles  were  affected  by  boundary  layers  near 
the  walls.  Nevertheless  it  would  appear  that  the  other  differences  observed 
between  his  experimental  and  theoretical  velocity  profiles  have  to  be  attributed 
to  the  fact  that  the  theoretical  result  is  based  on  a  periodic  solution. 

Faced  with  these  uncertainties,  many  designers  have  merely  sketched 
plausible  shapes  by  eye  or  selected  a  contour  used  previously  elsewhere  and 
believed  to  be  acceptable.  While  many  of  these  contractions  have  proved  to  be 
satifactory,  the  ability  to  consider  a  new  design  in  terms  of  a  theoretical  model 
must  be  considered  an  advantage. 

Batchelor  and  Shaw  (1944)  were  among  the  first  to  apply  numerical  methods 
to  examine  the  flow  properties  of  an  axisymmetric  contraction.  The  laborious 
computations  for  the  current  LSWT  contraction  were  performed  by  hand.  As 
mentioned  in  the  introduction  the  design  is  characterised  by  an  initially  rapid  rate 
of  contraction  followed  by  a  more  gradual  variation  near  the  exit.  The  results  of 
their  computations  led  them  to  suspect  that  the  inlet  boundary  layer  of  the  LWST 
may  suffer  from  intermittent  boundary  layer  separation  causing  increased 
unsteadiness  of  the  flow  in  the  working  section. 

Morel  (1975)  carried  out  a  study  aimed  at  providing  practical  guidelines  for 
the  design  of  axisymmetric  contractions.  Contours  formed  by  two  polynomial 
arcs  joined  smoothly  together  at  an  inflection  point  were  examined.  Cubic 
polynomials  were  found  to  form  the  most  suitable  wall  shape  and  these  were 
selected  for  a  detailed  study.  Finite-difference  solutions  of  the  Euler  equations 
were  used  to  produce  design  charts  which  define  the  contraction  length  and 
position  of  the  inflection  point  in  terms  of  allowable  adverse  pressure  gradients  at 
the  inlet  and  the  exit  and  the  desired  exit  flow  uniformity.  Morel  was  the  first  to 
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show  the  remarkable  fact  that  for  fixed  requirements,  the  contraction  length 
decreases  as  the  contraction  ratio  increases. 

Mikhail  (1979)  proposed  a  method  for  the  optimum  design  of  a  wind  tunnel 
contraction  defined  to  be  the  shortest  one  that  avoids  boundary  layer  separation 
and  which  supplies  the  flow  to  the  working  section  with  a  specified  uniformity. 
Boundary  layer  calculations  were  made  using  a  computer  program  based  on  the 
’lag-entrainment’  method  which  takes  into  account  the  wall  curvature.  He 
showed  that  by  optimising  the  contraction  wall  curvature  it  is  possible  to  reduce 
the  contraction  length  considerably.  Turbulent  flow  separation  was  predicted 
using  Stratford’s  criterion.  The  behaviour  of  turbulent  boundary  layers  with  extra 
strain  rates  (such  as  those  caused  by  streamline  curvature,  streamline  divergence 
or  pres?"re  gradients)  is  still  the  subject  of  current  research.  It  is  doubtful 
whether  Stratford’s  criterion  can  be  applied  with  any  confidence  to  boundary 
layers  on  a  contraction  wall  which  experience  streamline  curvature,  streamline 
convergence  and  pressure  gradients.  It  is  unfortunate  that  Mikhail  has  provided 
no  experimental  verification  for  his  suggestions.  Therefore  the  behaviour 
predicted  by  his  study  must  be  treated  with  caution. 

Recently  Downie,  Jordinson  and  Barnes  (1984)  described  a  method  for 
calculating  incompressible  inviscid  flow  through  ducts  of  rectangular  cross- 
section.  The  method  is  based  on  a  finite  difference  numerical  solution  of 
Laplace’s  equation  for  the  velocity  potential  in  three  dimensions.  Downie  et  al 
compared  the  prediction  of  their  numerical  calculations  with  experimental 
measurements  and  found  satisfactory  agreement.  The  magnitudes  of  the  pressure 
overshoots  and  undershoots  were  found  to  be  smaller  in  the  experiments 
presumably  because  of  the  cushioning  effects  of  boundary  layer  thickening.  An 
interesting  aspect  of  their  work  is  the  plotting  of  contours  of  the  pressure 
coefficient  distribution  superimposed  on  the  sides  of  the  contraction  walls  (see 
their  figure  6).  These  results  have  certain  implications  regarding  the 
approximation  of  axial  symmetry  as  applied  to  ducts  of  octagonal  cross-section. 

For  ducts  of  rectangular  cross-section  they  found  that  the  regions  of  pressure 
overshoot  and  undershoot  are  characterised  by  closed  contours  located  near  the 
contractions  entry  and  exit  planes  and  which  are  centred  on  the  duct  comers.  It 
is  likely  that  in  an  octagonal-sec tioned  contraction  similar  overshoots  and 
undershoots  of  the  pressure  coefficient  would  occur  at  the  corners  and  these 

would  be  larger  than  those  seen  along  the  wall  centrelines.  Since  an  octagonal 
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section  is  a  better  approximation  to  a  circular  section,  than  is  a  rectangular 
section,  and  since  the  effects  of  boundary  layers  in  a  real  cont- action  appear  to 
reduce  the  pressure ,  overshoots  and  undershoots,  it  seems  reasonable  to 
approximate  an  octagonal  section  with  axial  symmetry. 

4.  EXPERIMENTAL  VERIFICATION  OF  NUMERICAL  METHOD 

A  numerical  method  has  been  developed  for  the  solution  for  Stokes’s  stream 
function  (i.e.  equation  i)  with  the  prescribed  boundary  conditions.  Since  the 
equation  is  elliptic  the  solution  is  affected  by  upstream  and  downstream 
conditions  so  the  computational  doman  is  fitted  with  parallel-walled  extensions  to 
the  contraction  inlet  and  outlet  that  are  both  0.225  inlet  diameters  long.  The 
numerical  techniques  are  described  in  detail  in  Appendix  A.  Values  of  the 
streamwise  velocity  U  and  the  radial  velocity  V  are  obtained  by  differentiating 
cubic  splines  that  are  fitted  to  the  solutions  tor  the  stream  function  at  each  grid 
point  and  to  the  boundary  values.  Pressure  coefficients  Cp  based  on  the 
streamwise  velocity  at  the  most  downstream  point  on  the  axis  are  calculated  from 
the  Bemouilli  equation  so  that  the  variation  of  the  Cp  along  the  contraction  walls 
may  be  examined.  Contours  of  constant  Cp  can  also  be  calculated. 

To  test  the  approximations  of  axial  symmetry  made  for  octagonal  sections 
and  for  the  effects  of  neglecting  viscous  forces,  the  pressure  coefficient 
distribution  for  the  existing  LSWT  contraction  shape  has  been  calculated  for 
comparison  with  exp  rimental  measurements.  The  equivalent  radius  of  the 
axisym metric  contraction  was  determined  from  the  octagonal  cross-sectional  area 
at  12  axial  positions.  The  dimensions  of  the  octagonal  sections  were  taken  from  a 
drawing  which  was  used  to  construct  a  1/3.44  scale  model  of  the  LSWT  known  as 
the  Pilot  tunnel.  The  calculated  wall  pressure  coefficients,  shown  in  figure  2,  are 
characterised  by  an  almost  constant  value  of  Cp  from  the  beginning  of  the 
parallel-walled  entry  section  (Cp  -  0.9375)  through  to  near  the  entry  of  the 
contraction  where  an  abrupt  change  to  a  region  of  almost  constant  pressure 
gradient  occurs  for  about  1/8  of  the  contraction  length. 

The  pressure  co-efficient  levels  off  to  a  maximum  value  of  0.966  before  falling  to 
a  minimum  value  of  -0.062  near  the  contraction  exit.  A  short  region  of  adverse 
pressure  gradient  follows  the  Cp  minimum  which  is  about  twice  the  magnitude  of 
the  adverse  pressure  gradient  near  the  inlet.  The  contours  shown  in  figure  3 
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indicate  that  the  influence  of  the  contraction  extends  further  upstream  than 
downstream.  At  the  end  of  the  contraction  the  entry  flow  to  the  working  section 
is  uniform  to  within  around  ±1%. 

Experimental  measurements  of  the  wall  Cp  distribution  have  been  made  in 
the  Pilot  tunnel  which  differs  from  the  LSWT  in  that  it  is  of  the  open  return 
type.  The  maximum  working  section  velocity  is  around  20  m/s.  Figure  4  is  a 
sketch  of  the  tunnel  showing  relevant  dimensions.  Unfortunately  the  entry  flow 
is  characterised  ly  an  unsteady  vortex  originating  on  the  floor.  An  effort  was 
made  to  reduce  the  strength  of  the  vortex  by  raising  the  tunnel  by  300mm  and  by 
moving  the  entrance  away  from  a  nearby  wall.  Another  cause  of  flow  problems  is 
that  the  tunnel  entrance  and  exit  are  located  in  different  rooms  which  are 
connected  by  an  opening  through  which  the  tunnel  passes.  The  return  flow  to  .he 
entrance  was  made  more  uniform  by  doubling  the  area  of  the  opening  and  by 
placing  baffles  at  locations  selected  by  trial  and  error.  Observations  of  smoke 
released  on  the  floor  near  the  entrance  indicated  that  the  strength  of  the  vortex 
had  been  reduced  considerably.  A  weak  vortex  remained  which  appeared  and 
disappeared  sporadically. 

Static  pressure  tappings  were  placed  along  the  contraction  walls  at  equal 
axial  intervals.  Five  rows  of  46  tappings  each  are  located  along  Uie  centrelines 
of  the  contraction  top  wall,  fillet  and  side  wall  and  along  the  two  vertices 
between  these  three  surfaces,  within  one  quadrant  of  the  octagonal  section  as 
shown  in  figure  5.  All  tappings  are  connected  to  48  port  solenoid-driven 
scanivalves  which  are  driven  under  computer  control.  A  Pitot-static  tube  is 
fitted  in  the  working  section  to  provide  the  reference  static  and  total  pressures. 
The  results  of  the  experimental  Cp  measurements  are  shown  in  figures  6(a)  to 
(e).  All  distributions  have  the  characteristic  Cp  overshoot  near  the  contraction 
entrance  and  an  undershoot  near  the  exit.  Slightly  larger  undershoots  are 
observed  for  the  comer  distributions  which  may  be  due  to  the  corner  vortices. 
The  distributions  on  the  top  and  side  walls  and  along  comer  A  show  small  positive 
excursions  after  the  expected  undershoot  near  the  contraction  exit.  These 
excursions  may  be  due  to  small  misalignments  at  the  junction  between  the 
contraction  and  working  section. 

All  the  experimental  measurements  and  the  results  of  the  numerical 
calculation  are  superimposed  in  figure  7.  The  agreement  of  the  experimental 
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and  numerical  results  is  reasonable  except  for  the  distribution  measured  along  the 
top  wall.  The  top  wall  Cp  distribution  has  the  lowest  values  in  the  first  half  of 
the  contraction  while  it  has  the  highest  values  in  the  second  half,  yet  the  entrance 
and  exit  values  are  in  agreement  with  the  other  distributions.  These  results  may 
be  due  to  large-scale  geometrical  inaccuracies  in  the  contraction  shape  here  or  to 
the  effect  of  the  entry  vortex.  The  distributions  measured  on  the  -..ntrelines  of 
the  top,  fillet  and  side  walls  have  smaller  Cp  undershoots  than  that  predicted  by 
the  numerical  calculation.  This  is  consistent  with  the  results  of  Downie  et  al  who 
found  that  the  experimental  overshoots  and  undershoots  were  less  than  those 
predicted  presumably  because  of  local  boundary  layer  thickening.  On  the  basis  of 
the  results  presented  here  it  appears  that  the  approximation  or  circular  symmetry 
for  an  octagonal-sectioned  contraction  is  reasonable  and  that  the  numerical 
calculations  for  a  new  LSWT  contraction  shape  will  act  as  a  satisfactory  guide  for 
predicting  its  performance. 

5.  SELECTION  OF  A  NEW  CONTRACTION  SHAPE  FOR  THE  LSWT 


A  reasonable  criterion  for  the  selection  of  a  new  contraction  shape  for  the 
LSWT  is  to  minimise  the  adverse  pressure  gradient  near  the  inlet  while  not 
significantly  increasing  the  adverse  pressure  gradient  near  the  exit.  A  variety  of 
contraction  shapes  can  be  specified  by  two  simple  polynomial  arcs  which  are 
joined  together  at  an  inflection  point  (see  Morel  1975).  Using  the  notation  of 
fig-  ~e  8  then, 

for  0<x<Xj,  R  =  R^-AjXn  2(a) 

and  for,  X1<x<L(;,  R  *  Ro  +  A2  [Lc  -  x]n  2(b) 


where  n  is  an  exponent  which  controls  the  maximum  steepness  of  the  contraction 
shape.  By  matching  the  gradient  of  the  arcs  at  the  inflection  point  the  constants 
A^  and  Ag  are  given  by, 
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From  intuition  one  would  expect  that  the  position  and  magnitude  of  the 
adverse  pressure  gradients  near  the  inlet  and  outlet  of  a  contraction  to  be  related 
to  the  wall  curvature  distribution.  Morel  tested  six  power-law  contraction  curves 
of  order  n  =  2,  2.5,  3.0,  3.5,  and  4.0,  all  of  which  had  the  same  contraction  ratio 
(9:1),  length  to  inlet  diameter  ratio  (l_c/D^  =  1)  and  nondimensional  distance 
from  the  inlet  to  the  inflection  point  (X=X^/LC=0.5) .  As  the  exponent  n  is 
reduced  the  wall  curvature  near  both  the  entry  and  the  exit  decreases  so  that  the 
pressure  extrema  and  adverse  pressure  gradients  also  may  be  expected  to 
decrease.  However  Morel’s  computations  indicate  that  this  effect  only  occurs  for 
exponents  greater  than  3. 

To  check  that  this  surprising  result  also  occurs  for  the  contraction  ratio  of 
the  LSWT  the  Cp  distributions  for  the  contraction  curves  given  by  n  =  2,  3,  4  and  5 
were  calculated.  All  the  contraction  curves  had  the  same  nondimensional 
distance  to  the  inflection  point  (X  =  0.6)  and  length-to-inlet-diameter  ratio  of  the 
LSWT  ( Lc/D .  =  0.881).  The  wall  Cp  distributions  shown  in  figure  9  indicate 
that  as  the  polynomial  order  is  increased  from  3  to  5  both  the  entry  and  exit  Cp 
extrema  and  adverse  pressure  gradients  increase.  Although  the  Cp  maximum  is 
small  for  n  =  2,  the  adverse  pressure  gradient  is  still  relatively  large.  For  n  =  2 
the  Cp  minimum  near  the  exit  is  larger  than  for  n  =  3  and  the  adverse  pressure 
gradient  is  larger  than  all  the  other  contraction  curves.  The  ratios  of  the  adverse 
pressure  gradients  to  those  estimated  from  the  LSWT  contraction  calculation  are 
approximately  0.40,  0.25,  0.36  and  0.67  near  the  entry  and  4.3,  2.3,  2.6  and  3.3 
near  the  exit  for  the  contraction  shapes  given  by  n  =  2,  3,  4  and  5  respectively. 
The  Cp  contours  for  n  =  2,  3  and  4  shown  in  figures  10(a)  to  (c)  indicate  that  the 
upstream  and  downstream  influences  of  the  contraction  become  larger  as  n  is 
reduced.  In  all  cases  the  greatest  proportion  of  the  working  section  length  is 
required  on  the  axis  for  the  velocities  to  rise  to  within  about  1%  of  the  downsteam 
reference  velocity.  Less  working  section  length  is  required  for  the  velocity 
overshoot  at  the  boundary  to  fall  to  within  the  same  tolerance. 

Compared  to  the  existing  LSWT  contraction  the  third-order  contraction 
shape  has  an  inlet  adverse  pressure  gradient  which  is  smaller  by  a  factor  of  4 
while  the  adverse  pressure  gradient  near  the  exit  is  larger  by  a  factor  of  2.3. 
Also  the  favourable  pressure  gradient  near  the  centre  of  the  contraction  is  about 
twice  that  of  the  LSWT  so  that  there  is  the  possibility  of  boundary  layer 
relaminarisation  here.  Since  laminar  layers  are  more  prone  to  separate  than  are 
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turbulent  layers,  the  increased  adverse  pressure  gradient  at  the  exit  may 
introduce  new  problems.  It  appears  that  the  cubic  profile  with  X  =  0.6  may 
produce  too  large  a  modification  to  the  adverse  pressure  gradients.  If  the 
exponent  n  in  equations  2(a)  and  (b)  is  held  constant  the  effect  of  moving  the 
inflection  point  towards  the  contraction  entry  would  be  to  decrease  the  curvature 
near  the  exit  and  increase  the  curvature  near  the  entry.  Hence  the  Cp  minimum 
and  the  adverse  pressure  gradient  near  the  exit  would  decrease  at  the  expense  of 
increasing  the  Cp  maximum  and  adverse  pressure  gradient  near  the  entry. 
Therefore  solutions  were  obtained  for  the  cubic  contraction  shape  with  X  =  0.5 
and  l_c/D..=0.881 .  Calculations  were  also  performed  for  contraction  shapes 
with  n  =  2.5  and  n  =  2.75  to  further  examine  the  curious  affect  noted  earlier 
regarding  the  Cp  near  the  exit  for  n  <  3. 

The  wall  Cp  distributions  are  shown  in  figure  11.  The  adverse  pressure 
gradients  near  the  contraction  exit  are  all  of  similar  magnitude.  However,  the 
contraction  shape  given  by  n  =  3  offers  the  smallest  Cp  undershoot  and  adverse 
pressure  gradient  near  the  exit.  The  contours  for  the  case  n  =  3  are  shown  in 
figure  12  and  indicate  that  the  entry  flow  to  the  working  section  approaches 
uniformity  more  rapidly  than  the  case  where  n  =  3  and  X  =  0.6  (see  figure  10b). 
As  may  be  expected  from  the  wall  curvature  distribution,  the  upstream  influence 
of  the  contraction  is  greater  for  the  case  X  =  0.5  than  for  X  =  0.6.  The  entry/exit 
pressure  gradients  are  approximately  0.38  and  1.65  times  those  calculated  for  the 
LSWT  contraction. 

It  appears  that  to  achieve  the  criterion  stated  earlier  for  a  new  LSWT 
contraction  it  is  necessary  to  increase  the  contraction  length.  Modem 
contraction  designs  usually  have  a  length  to  inlet  diameter  ratio  around  1  (see 
Mikhail  1979).  A  longer  contraction  could  easily  be  accomodated  by  the  LSWT 
since  it  has  a  rather  longer  settling  chamber.  The  calculated  wall  Cp 
distribution  for  the  contraction  shape  specified  by  n  =  3,  X  =  0.5  and  Lc/D.=  1  is 
shown  in  figure  13.  The  inlet  and  outlet  pressure  gradients  are  approximately 
0.20  and  0.91  times  those  calculated  for  the  LSWT  contraction.  This  contraction 
satisfies  the  criterion  of  reducing  the  entry  adverse  pressure  gradient  (by  a  factor 
of  5)  while  not  increasing  the  adverse  pressure  gradient  near  the  exit.  The 
pressure  coefficient  contours  shown  in  figure  14  indicate  that  the  upstream  and 
downstream  influence  of  this  contraction  shape  is  small.  In  particular  the  flow  at 
the  exit  of  the  contraction  is  uniform  to  within  around  ±1%. 
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6.  CONCLUDING  REMARKS 

To  avoid  interference  with  the  steel  shell  the  new  design  would  have  to  be 
fitted  into  the  LSWT  upstream  of  the  existing  contraction.  The  new  contraction 
could  be  made  in  the  form  of  a  fibreglass  covered  plywood  skin  which  is  supported 
by  timber  beams  with  much  of  the  material  being  pre-fabricated  so  that  the 
installation  would  proceed  quickly.  The  results  of  the  calculations  suggest  that  a 
much  shorter  settling  section  could  be  used  without  affecting  the  uniformity  of 
the  entry  flow  to  the  contraction  (see  figure  A2).  According  to  the  computations 
the  settling  region  could  be  reduced  in  length  by  as  much  as  10ft.  (*3m> 

leaving  room  for  a  screen  to  be  installed  after  the  honeycomb.  The  suggested 
arrangement  shown  in  figure  15  would  provide  an  extra  8  ft.  (»2.4m)  of  working 
section. 
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APPENDIX  A  NUMERICAL  TECHNIQUES 

A.l  Relaxation  method  for  solution  of  Stokes’s  stream  function. 

A  calculation  method  has  been  developed  by  author  for  the  solution  of 
Stokes’s  Stream  Function,  i.e.  equation  (1).  The  code  was  written  in  FORTRAN 
77  and  was  developed  and  run  on  a  16-bit  PDP  11  /23-PLUS  microcomputer.  The 
development  of  the  code  on  the  microcomputer  allowed  rapid  turnaround  time 
during  the  edit,  compile,  task-build  and  debug  stages  of  program  development 
without  incurring  the  inconvenience  of  mainframe  computing.  One  of  the 
problems  of  using  a  16-bit  computer  for  numerical  computations  is  the  small 
maximum  program  size  available  (65,536  16-bit  words)  which  is  much  less  than  the 
typical  memory  size  available  with  these  machines.  In  order  to  make  use  of  extra 
memory  PDP  11  FORTRAN  77  allows  arrays  to  be  declared  as  VIRTUAL  meaning 
that  array  values  are  stored  outside  the  65K  program  boundary.  Hence  large 
arrays  do  not  consume  program  space  which  can  be  used  for  extra  computational 
code. 


Integration  of  the  governing  differential  equation  with  the  prescribed 
boundary  conditions  is  accomplished  using  a  numerical  method  based  on  the 
calculus  of  finite  differences  known  as  the  relaxation  method.  Details  of  the 
method  can  be  found  in  undergraduate  textbooks  (e.g.  Smith  1978).  The 
computational  grid  is  generated  by  dividing  the  rectangular  region  up  into  equal 
sized  intervals  in  the  axial  and  in  the  radial  directions.  Currently  the  maximum 
number  of  grid  lines  in  either  direction  is  one  hundred  giving  a  maximum  number 
of  nodal  intersection  points  of  10,000.  For  all  cases  reported  here  equal  numbers 
of  grid  lines  are  used  in  both  the  axial  and  radial  directions.  Since  the  axial 
length  of  the  computational  region  is  always  larger  than  the  radial  height  much 
finer  grid  resolution  is  obtained  in  the  radial  direction.  The  finer  resolution  in 
the  radial  direction  is  an  advantage  due  to  the  larger  gradients  of  the  stream 
function  experienced  near  the  tunnel  boundaries.  No  attempt  was  made  to  use  a 
non-uniform  grid  to  cope  with  the  larger  gradients  of  the  stream  function 
experienced  near  the  tunnel  walls,  nor  was  a  coordinate  transformation  utilised  to 
make  the  calculation  domain  rectangular  (e.g.  Mikhail  1979)  since  the  ratio  of  the 
inlet  to  outlet  diameters  is  only  2:1.  Initial  values  are  assigned  to  the  nodal 
points  as, 

*  =  <r/rw>2 
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These  values  are  systematically  adjusted  to  satisfy  equation  (1)  and  the  boundary 
conditions  in  an  interative  manner.  Nodal  points  outside  the  tunnel  boundary  are 
not  used  in  the  computation.  Figure  A1  shows  a  grid  region  where  the 
contraction  boundary  intersects  the  grid  points  at  A  and  B.  In  this  region  the 
lengths,  a,  b,  c  and  d  are  not  equal.  Referring  to  figure  A1  it  can  be  Bhown  that, 
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Further  within  the  contraction  at  points  A  ,  B’,  C’,  and  D’  for  example,  the  lengths 
a’  and  c’  are  equal,  and  b’  and  d’  are  equal  so  that  the  expression  for  ip  above  is 
simplified. 

Each  iteration  consists  of  calculating  the  new  estimate  (<PQ)  for  each  nodal 
point,  except  those  in  the  first  and  last  columns,  using  the  values  from  the 
preceding  iteration.  The  values  for  the  first  column  specify  the  uniform  and 
parallel  entry  flow  so  that  the  values  here  are  kept  fixed.  The  values  for  the 

points  on  the  last  column  are  evaluated  from  the  finite  difference  approximation 

2  2 

to  the  boundary  condition  3  ifi/dx  -0.  The  node  with  the  maximum  percentage 
change  from  the  previously  iterated  value  is  stored  during  each  new  iteration 
cycle.  The  iteration  process  is  halted  when  this  maximum  percentage  change 
falls  below  a  user  specified  value.  For  all  the  results  presented  here  the 
iterations  were  continued  until  all  nodal  values  remained  unchanged  from  their 
previous  value  i.e.  the  results  had  converged  to  around  7  significant  figures. 
Although  less  convergence  could  be  tolerated  the  author  did  not  wish  to  become 
involved  with  deciding  what  a  lesser  tolerable  convergence  level  might  be. 

The  input  to  the  program  is  an  ASCII  data  file  containing  pairs  of  values 
representing  the  axial  and  radial  coordinates  of  the  tunnel  profile  to  be  tested. 
The  origin  of  these  data  pairs  (whether  scaled  from  a  sketch  or  obtained  from  a 
specified  function)  does  not  matter  so  long  as  the  radius  is  constant  or  decreasing 
with  streamwise  distance.  It  is  the  user’s  responsibility  to  include  adequate 
entry/exit  region  lengths  in  the  file.  A  cubic  spline  is  fitted  to  these  data 
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points.  The  grid  is  generated  on  a  rectangular  region  of  length  equal  to  the  axial 
distance  between  the  first  and  last  data  pairs  and  of  height  1.001  times  the  inlet 
radius.  The  multiplying  factor  (1.001)  is  used  to  ensure  that  the  tunnel  boundary 
falls  wholly  within  the  rectangular  region  and  that  the  tunnel  boundary  values 
would  fall  slightly  over  and  some  slightly  under  a  grid  line  causing  computational 
difficulties.  It  was  found  that  the  total  computational  time  varied  approximately 
as  (where  N  =  the  number  of  axial  grid  lines  =  the  number  of  radial  grid  lines). 
After  some  initial  trials  it  was  found  that  N  =  60  would  give  results  of  ample 
resolution.  For  N  =  60  the  computations  take  around  15  hours  (i.e.  an  overnight 
run)  performing  around  4000  iterations.  At  the  end  of  each  run  a  data  file  is 
created  for  subsequent  analysis. 

A. 2  Evaluation  of  wall  pressure  coefficients. 

The  steam  function  solutions  are  processed  to  produce  a  file  of  the  wall 
pressure  coefficients  for  each  axial  grid  coordinate.  A  cubic  spline  is  applied  to 
the  stream  function  solutions  at  each  grid  point  and  at  the  boundary  along  each 
column  in  succession.  Values  of  U  (=  1/rdp/dr)  are  found  at  the  boundary  from 
the  derivation  of  the  spline  fit.  Since  the  boundary  is  a  streamline  and  its  slope 
can  be  calculated  from  the  splined  values  of  the  boundary  shape,  the  radial 
velocity  component  can  also  be  calculated  at  the  boundary.  The  values  of  the 
pressure  coefficient  at  each  columnar  grid  boundary  position  are  simply  evaluated 
from  the  Bernoulli  equation  as, 

Cp  =  1  -  <UB2  +  VB2>  /  U  REF2  <A» 

where  Ug  and  Vg  are  the  streamwise  and  radial  boundary  velocity  components  and 
UreF  ’8  t*le  streamwise  velocity  of  the  most  downstream  axis  grid  point.  Values 
of  U  are  indeterminate  on  the  axis  (since  U  =  1/r  dip /dr  =  0/0).  A  fifth  order 
least  squares  polynomial  of  best  fit  passing  through  the  closest  point  to  the  axis  is 
extrapolated  to  the  axis  to  evaluate  Uggp. 
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A.3  Contour  plotting. 

Contours  of  constant  pressure  coefficient  are  calculated  using  a  second- 
order  predictor-corrector  scheme  which  requires  evaluation  of  the  stream  wise  and 
radial  directional  derivatives  of  surfaces  describing  the  Cp  distribution  over  each 
grid  region.  The  Cp  values  are  determined  from  the  Bemouilli  equation  (using 
equation  A2  with  the  same  reference  velocity  but  using  each  grid  value  of  U  and  V 
instead  of  the  boundary  values).  Internal  values  of  U  and  columnar  boundary 
values  of  U  and  V  are  determined  using  the  same  techniques  as  described 
earlier.  Similarly,  cubic  splines  are  applied  to  the  stream  function  along  each 
row  (excluding  the  axis)  for  internal  values  of  V  at  each  grid  point  and  row 
boundary  values  of  V  and  U  (using  the  slope  of  the  boundary).  Axis  values  of  V 
are  zero  by  definition.  Therefore  values  of  U  and  V  (and  hence  Cp)  can  be  found 
at  each  grid  point  inside  the  contraction  boundary  and  at  the  intersection  of  each 
column  and  row  with  the  boundary. 

Over  each  grid  region  through  which  the  contour  is  plotted,  a  surface 
polynomial  in  x  and  r  is  used  to  describe  the  variation  of  Cp  It  is  essential  for 
accurate  contouring  that  each  surface  possess  common  boundaries  with  those  of 
the  four  adjoining  regions.  These  boundaries  are  the  result  of  cubic  splines 
applied  along  the  appropriate  rows  and  columns.  Since  cubic  splines  are  used 
each  of  the  common  boundaries  can  be  fitted  exactly  with  a  third  order 
polynomial  By  adding  two  fourth  order  terms,  linear  in  x  and  r  respectively,  the 
requirements  of  the  surface  fit  can  be  satisfied  i.e. 


+  AgX  +  A^r 

+  A.x2  +  A,x  r  +  A,r2  (A3) 

4  5  6 

+  A;x3  +  A8x2r  +  Agx  r2  +  A10r3 
+  Allx  r3  +  A12x  r3 

The  12  unknown  coefficients  A^  ...  Aj^  are  evaluated  from  12  equation*  formed 
from  the  following  conditions  :  the  surface  must  pass  through  the  four  vertices  of 
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the  grid  (4  equations)  and  the  first  derivatives  aCp/ax  and  3Cp/ar  at  each  of  the 
vertices  must  be  equal  to  that  Indicated  by  the  cubic  splines  <8  equations).  While 
there  is  continuity  of  the  Cp  values  across  each  boundary,  continuity  of  first 
derivatives  cannot  be  guaranteed  except  at  the  comers.  The  12  coefficients  can 
be  obtained  by  the  inversion  of  a  12x12  matrix.  When  the  contour  moves  into  a 
new  region  a  new  set  of  12  coefficients  must  be  used.  In  order  to  save 
computation  time  a  localised  stretched  coordinate  system  x1  r1  is  used  the 
origin  at  the  centre  of  the  grid  region  so  that  the  region  boundaries  line  on  the 
lines  x'  =  ±1,  r'  =  +1.  The  results  of  the  inversion  of  the  12x12  matrix  below 
can  be  performed  once  and  stored  so  that  the  contouring  program  can  rapidly 
calculate  the  coefficients  A1  ...  A12  for  each  new  grid  region  as  required. 
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Problems  arise  with  the  surface  fitting  technique  since  grid  regions  which 
are  cut  by  the  tunnel  boundary  require  values  and  derivatives  of  Cp  that  are 
outside  the  boundary  i.e.  some  form  of  extrapolation  is  required.  Depending  on 
the  relative  size  of  the  axial  and  radial  grid  intervals  and  on  the  slope  of  the 
contraction  wall,  regions  were  found  where  points  up  to  three  radial  grid  intervals 
outside  the  boundary  were  required  for  the  surface  fitting.  Fifth-order  least 
squares  polynomials  passing  through  the  boundary  values  ire  used  to  the 
extrapolation.  Purely  radial  extrapola  ons  are  used  where  the  boundary  radius  is 
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greater  than  0.9  times  the  inlet  radius  or  less  than  1.1  times  the  exit  radius.  For 
extrapolation  where  the  boundary  radius  falls  within  these  limits,  the  points 
outside  the  boundary  that  directly  influence  the  contours  are  calculated  using  both 
radial  and  streamwise  extrapolations.  The  extrapolated  values  are  averaged  by 
weighting  them  with  the  distances  from  the  grid  point  to  the  opposite  boundary 
respectively.  Fifth-order  polynomials  in  the  radial  direction  that  pass  through 
these  weighted  averaged  points  are  used  to  extrapolate  for  the  values  of  U  and  V 
at  the  remaining  grid  points.  Values  are  required  at  these  points  so  that  the  grid 
region  boundary  splines  in  the  regions  used  for  contouring  are  not  adversely 
affected. 

A. 4  Selection  of  inlet  and  outlet  extension  lengths. 

Since  the  governing  equation  is  elliptic  in  nature  the  flow  within  the 
contraction  is  sensitive  to  the  upstream  and  downstream  flow  conditions. 
Therefore  it  is  important  that  the  lengths  of  the  straight-walled  sections  before 
and  after  the  contraction  are  large  enough  that  truncation  beyond  them  has 
negligible  influence  on  the  results.  Morel  (1976)  used  inlet  and  outlet  extensions 
that  were  0.7  local  diameters  in  length.  He  reports  that  further  lengthening  had 
no  influence  on  his  results  but  made  no  mention  of  the  effects  of  shortening  the 
extensions.  A  series  of  calculations  have  been  made  using  equal  inlet/outlet 
extensions  of  lengths  0.125,  0.225,  0.325  and  0.425  times  the  inlet  diameter  for 
the  cubic  contraction  shape  shown  in  figure  10(b).  The  Cp  contours  shown  in 
figures  A2(a)  to  (d)  reveal  that  only  the  contours  near  the  inlet  and  outlet  of  the 
contraction  with  the  nondimensional  inlet/outlet  extension  length  of  0.125  are 
significantly  affected.  Differences  between  the  contours  for  the  other  cases  are 
undetectable.  Therefore  inlet/outlet  extensions  that  are  0.225  inlet  diameters 
long  were  used  for  all  computations  unless  otherwise  specified. 

APPENDIX  B  EXPERIMENTAL  MEASUREMENTS  OF  PRESSURE  COEFFICIENT. 

Pressure  coefficients  have  been  measured  along  the  contraction  walls  of  the  Pilot 
(LSWT  model)  tunnel.  Five  rows  of  46  static  pressure  tappings  each  are  located 
along  the  centrelines  of  the  contraction  top  wall,  fillet  and  side  walls  and  along 
the  two  vertices  between  these  three  surfaces,  within  one  quadrant  of  the 
octagonal  section  as  shown  in  figure  5.  The  tappings  are  positioned  according  to 
arc  lengths  which  has  been  calculated  from  the  contraction  geometry.  The 
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calculated  arc  lengths  over  the  contraction  length  (around  1500mm)  were  found  to 
be  within  ±lmm  of  the  measured  values.  Therefore  the  axial  positions  of  the 
tappings  can  be  regarded  as  precise.  All  tappings  are  connected  to  48  port 
solenoid-driven  scanivalves  via  plastic  tubing.  Only  three  scanivalves  are  used 
for  simultaneous  measurements  throughout  the  experiments  due  to  a  lack  of 
pressure  transducers  with  a  suitable  range.  The  three  external  differential 
pressure  transducers  used  are  : 

(1)  Setra  model  236  with  a  range  of  0.2  psi  and  full  scale  output  of  5v. 

(2)  Setra  model  236  with  a  range  of  0.5  psi  and  full  scale  output  of  5v. 

(3)  Furness  FC  6124  with  a  range  of  1.5  psi  and  full  scale  output  of  lv. 

The  outputs  of  the  transducers  are  suitably  amplified  before  measurement.  The 
two  extra  ports  on  each  scanivalve  are  connected  to  total  and  static  pressures  of  a 
Pitot-Static  tube  located  in  the  working  section.  The  static  pressure  from  the 
Pitot-Static  tube  is  also  used  as  a  reference  pressure  for  the  transducers. 
Transducer  drift  is  monitored  and  corrected  for  during  the  experiments  by 
connecting  the  reference  pressure  across  each  transducer.  Since  five  groups  of 
pressure  tappings  are  measured,  each  experimental  run  was  performed  by 
connecting  sets  of  three  and  two  groups  of  connectors  to  the  scanivalves  in 
succession. 

The  scanivalve  solenoids  are  operated  by  a  controller  with  two  inputs, 
namely  ’single  step’  and  ’return  to  home.’  The  PDP  11/23-PLUS  microcomputer  is 
used  to  operate  the  controllers  and  to  store  the  amplified  transducer  outputs  via  a 
data  logger  with  the  brand  name  of  Dataporte.  The  Dataporte  is  a  compact 
intelligent  microprocessor-based  data  acquisition  system  which  can  be  supervised 
by  any  computer  with  an  RS232C  serial  interface.  Amongst  the  features  of  the 
Dataporte  are  analog,  digital,  event  and  counter  inputs,  data  storage  and 
averaging,  analog  and  digital  outputs,  temperature  measurement  and  a  real  time 
clock.  Up  to  46  single-ended  analog  input  channels  can  be  connected  with 
autoranging  over  three  decades.  The  Dataporte  is  controlled  by  issuing  simple 
abbreviated  English-like  commands.  The  maximum  sampling  rate  is  around  20 
samples/sec  so  that  the  system  is  only  suitable  in  applications  requiring  low 
sampling  rates.  One  of  the  eight  digital  outputs  is  used  to  step  the  three 
scanivalve  solenoids  in  parallel  at  a  rate  of  about  3  steps/sec  while  another  digital 

output  is  used  to  simultaneously  step  the  scanivalves  to  the  home  position. 
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Electronic  integration  of  the  amplified  pressure  transducers  can  be  used  to  obtain 
true  mean  values  over  time  periods  of  up  to  three  minutes  in  duration.  For  the 
results  presented  in  this  report  an  integration  time  period  of  15  seconds  was 
used.  Two  digital  outputs  of  the  Dataporte  are  used  to  start  and  reset  the 
integrators  while  one  of  the  digital  inputs  is  used  to  detect  completion  of  the 
integration  cycle.  The  integrated  voltages  are  then  returned  to  the  computer  for 
processing  and  storage. 


FIG.  1  SKETCH  OF  THE  ARL  LOW  SPEED  WIND  TUNNEL 
(LSWT)  SHOWING  IMPORTANT  DIMENSIONS. 
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FIG.  2.  CALCULATED  WALL  Cp DISTRIBUTION 

FOR  THE  EXISTING  LSWT  CONTRACTION. 
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FIG.  3.  CALCULATED  CONTOURS  OF  CpWITHIN  THE  LSWT. 


SKETCH  OF  OPEN  RETURN  PILOT  TUNNEL  (LSWT  MODEL) 
SHOWING  RELEVANT  DIMENSIONS 
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FIG.  7.  CALCULATED  AND  EXPERIMENTAL  VALUES  OF  Cp.  SYMBOLS  AS  IN 
FIGS.  2  AND  6.  EXCEPT  FOR  THE  EXPERIMENTAL  TOP  WALL 
MEASUREMENTS,  THE  CALCULATED  AND  EXPERIMENTAL  VALUES  ARE 
IN  REASONABLE  AGREEMENT. 
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FIG.  8.  NOTATION  FOR  CONTRACTION  WALL  CONTOUR 
CONSTRUCTED  FROM  TWO  ARCS. 


FIG.  9.  CALCULATED  WALL  Cp DISTRIBUTIONS  FOR  X=0.6,  Lc/D, =0.881 

(SAME  AS  LSWT)  AND  n=2,  3,  4  AND  5.  ARROW  SHOWS  ORDER  OF 
INCREASING  n. 
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FIG.  12.  CONTOURS  OF  Cp  WITHIN  CONTRACTION  SPECIFIED  IN  FIG.  11  WITH  n=3 
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FIG.  13.  CALCULATED  WALL  Cp DISTRIBUTION  FOR 
X=0.5,  Lc/Dj=1  AND  n=3. 
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FIG.  14.  CONTOURS  OF  Cp  WITHIN  CONTRACTION  SPECIFIED  IN  FIG  13. 

UNIFORMITY  OF  THE  ENTRY  FLOW  TO  THE  WORKING  SECTION  IS 
SIMILAR  TO  LSWT  CONTRACTION  (SEE  FIG.  3). 
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FIG.  15.  SKETCH  OF  EXISTING  SETTLING  REGION,  CONTRACTION, 

AND  WORKING  SECTION.  SUGGESTED  NEW  ARRANGEMENT 
PROVIDING  EXTRA  8  ft.  (2.4m)  OF  WORKING  SECTION  ALSO 
SHOWN. 


FIG.  A2.  CALCULATEDCpCONTOURS  FOR  DIFFERENT  INLET/OUTLET  EXTENSIONS  Xe. 
(a)  Xe=0.125Dj  (b)  Xe=0.225Dj  (c)  Xe=0.325Dj  (d)Xe=0.425D, 

ONLY  THE  CASE  WITH  X=0.125  HAS  SIGNIFICANTLY  DIFFERENT 
CONTOURS. 
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